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Recent studies of Charmonium Decays at CLEO 
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Department of Physics and Astronomy, University of Rochester, Rochester, New Yorl< 14627, USA 



Recent results on Charmonium decays are reviewed which includes two-, three- and four-body decays of x^-j 
states, observations of Y(4260) through ttttJ/i/i transitions, precise measurements of M{D^), M{rj) as well as 
Biv-'X). 
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1. Introduction 

Decays of a bound state of a quark and its anti- 
quark, quarkoniuni, provide an excellent laboratory 
for studying QCD. Particularly, heavy quarkonia such 
as charmonium states are less relativistic, thus play a 
special role in probing strong interactions. 

CLEO recently has accumulated data taken at the 
'0(25') resonance, providing a total of 27M ip{2S) de- 
cays. With the combination of this large statistical 
sample and the excellent CLEO detector, we will ex- 
plore an unprecedented world of charmonia. While 
many analyses are currently being carried out, in this 
note we present recent results on multi-body XcJ de- 
cays which employed the pre-existing 3M iJj{2S) sam- 
ple. 

We also present recent studies on decays of one of 
the exotic states, Y(4260), as well as precision mea- 
surement on AI(D^) that has an implication on prop- 
erties of X(3872). 

Finally, based on the full sample of ip{2S) data, we 
have results on properties of one of the light mesons, 
r]. 



combinations of rj and rj' . No Xci is seen as expected 
from conservation of spin-parity. 



0990906-0023 




3.30 



3.40 



ri'^'ri^'' Mass 



3.50 



3.60 



o 



X 



2. Factory of Xci^^Pj) states 

Xci^^Pj) states, which have one unit of orbital an- 
gular momentum and total spin of J=0, 1, or 2, can- 
not be produced directly from e^e~ collisions. They 
can be reached from ip{2S) through radiative (elec- 
tric dipole) transitions. Since B{'il;{2S) -^ IXc.i) = 
(9.3 ± 0.4, 8.8 ± 0.4^nd 8.1 ± 0.4) x IQ-^ for J=0, 1, 
and 2 respectively [l|, 27M ijj{2S) decays of the new 
data provides ~2M decays of each spin state of XcJ 
which should give us a greater understanding of the 
decay mechanisms of the XcJ rnesons. 

In this section, we present recent results of studies 
of XcJ decays based on 3M i^{2S) decays which should 
serve as the foundation for the future precision mea- 
surements by employing the full data sample of 27M 
of iIj{2S) decays. 



2.1. Two-body decay 

We present results on Xcj decay into combinations 
of ?7 and rj' mesons. Figure [T] shows invariant masses of 



Figure 1: Invariant masses of tjtj (a), rj'ri (b), and g'rj' (c) 



We measured B{xcO -^ VV) to be (0.31 ± 0.05 ± 
0.04 ± 0.02)% [2| where the first uncertainty is sta- 
tistical, the second is systematic, and the third is 
systematic due to the uncertainty in B{'ip{2S) -^ 
IXc.])- This is slightly higher, but consistent with, 
the two previously published measurements. The BES 
Collaboration measured this branching ratio to be 
(0.194 ± 0.085 ± 0.059)% f\ and the E-835 Collab- 
oration ^ had (0.198 ± 0.068 ± 0.037)%. We also 
measured B{xcO ^ v'v') to be (0.17 ± 0.04 ± 0.02 ± 
0.01)% for the first time. We set upper limits for 
B{Xco -^ m') < 0.05% , Bixc2 ^ w) < 0.047%, 
B{xc2 -^ W) < 0.023%, and B{xc2 -> v'v') < 0.031% 
at 90% confidence level. 

Our result can be compared to predictions based on 
the model of Qiang Zhao . He translates these decay 
rates into a QCD parameter, r, which is the ratio 
of doubly- to singly-OZI suppressed decay diagrams. 
In his model, our results indicate that the singly-OZI 
suppressed diagram dominates in these decays. 
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2.2. Three-body decay 

We have also looked at three-body deca;^ of XcJ 
states (one neutral and 2 charged hadrons) q]. They 
are tt+tt"/;, K^K~ri, pprj, TT~^Tr~rj' , K^K~t:^ , Ppi^^ ■, 
TT^K^Kg, and K^pK. Measured branching fractions 
are summarized in Table[Tl Again, our results are con- 
sistent with the results from BES Collaboration [7|, 
with better precision. 

In three of the above modes we looked for, 7r~''7r~77, 
i^T+iiT ~7r", and ■k^K~ K^, we observed significant sig- 
nals of Xci decays which are shown in Figure [21 
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Figure 2: Invariant masses of 7r^7r~r7 (top-left), K'^ K^-k 
(top-riglit), and K''^K^it^ (bottom) 

We performed Dalitz plot analyses based on these 
3 modes in which we neglected any possible interfer- 
ence effects between resonances and polarization of 
Xci- We estimated there could be ~ 20(15)% varia- 
tions in fit fractions for the mrri (KKn) mode due to 
such a simplified model. Figure [3] shows the Dalitz 
plot for TT^TT~r] and Table HTl shows its resultant fit 
fractions for each source. It is interesting to note that 
our data demand a relatively large yield of a cr pole. 

As for the KKtt mode, we performed simultaneous 
fits between Xci -^ K~^K~tt^ and Xci — >■ KsKtt by 
taking advantage of isospin symmetry. Dalitz plots for 
these modes are shown in Figures |4] and [S] Table IIIII 
shows their resultant fit fractions. 



2.3. Four-body decay 

We present a preliminary result on four-body decay 
of XcJ states in which we reconstructed h^h~Tr^Tr^, 
where h — tt, K, p; K+K~r]'!T°] and K'^tt^K^w'^. 
Results of this kind of study, many-body decays of 
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Figure 3: Dalitz plots for xa ^^ tj-k^tt 
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Figure 4: Dalitz plots for Xci — * K K tv 



XcJ states, should help to build a comprehensive un- 
derstanding about the P-wave dynamics. 

Clean signals were seen in all modes except Xci —^ 
ppTr'^TT^ for the first time as can be seen in Figure [6] 
Many resonant substructures were also seen for which 
we only considered significant ones (Tr+Tr^Tr^Tr^ and 
K^n^K'^TT^). The results are summarized in Table IVl 



Proceedings of the CHARM 2007 Workshop, Ithaca, NY, August 5-8, 2007 



Table I Branching fractions in units of 10~ . Uncertainties are statistical, systematic due to detector effects plus 
analysis methods, and a separate systematic due to uncertainties in the ^(2S) branching fractions. Limits are at the 
90% confidence level. 



Mode 


XcO 


Xcl 


Xc2 


n ir~ri 


<0.2f 


5.0 ±0.3 ±0.4 ±0.3 


0.49 ±0.12 ±0.05 ±0.03 


K+K-ri 


<0.24 


0.34 ±0.10 ±0.03 ±0.02 


< 0.33 


PPV 


0.39 ±0.ff± 0.04 ±0.02 


<0.f6 


0.19 ±0.07 ±0.02 ±0.01 


TT+TT-V 


<0.38 


2.4 ±0.4 ±0.2 ±0.2 


0.51 ±0.18 ±0.05 ±0.03 


K+K-7v° 


<0.06 


1.95 ±0.16 ±0.18 ±0.14 


0.31 ±0.07 ±0.03 ±0.02 


PpTv" 


0.59 ±0.fO± 0.07 ±0.03 


0.12 ±0.05 ±0.01 ±0.01 


0.44 ±0.08 ±0.04 ±0.03 


n+K-K° 


<0.10 


8.1 ±0.6 ±0.6 ±0.5 


1.3 ±0.2 ±0.1 ±0.1 


K+pA 


f.07±0.f7±0.f0±0.06 


0.33 ±0.09 ±0.03 ±0.02 


0.85 ±0.14 ±0.08 ±0.06 



Table II Fit results for Xci — * r]ii^-K~ Dalitz plot analysis. 
The uncertainties are statistical and systematic. Allowing 
for interference among the resonances changes the fit frac- 
tions by as much as 20% in absolute terms as discussed in 
the text. 



Mode 


Fit Fraction (%) 


ao(980)±7rT 


75.1 ±3.5 ±4.3 


/2 (1270)77 


14.4 ±3.1 ±1.9 


ar) 


10.5 ±2.4 ±1.2 
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Figure 5: Dalitz plots for Xci -^ irKKg. 



The measured branching fraction of Xcj -^ p tt^tt^ 
is consistent with that of Xcj -^ p^tt+tt" as ex- 
pected from isospin symmetry. Similar isospin sym- 
metry is also seen in Table ITVl where the partial width 



Table III Results of the combined fits to the Xci —^ 
K'^ K~-k'^ and Xa — > nKKg Dalitz plots. Allowing for in- 
terference among the resonances changes the fit fractions 
by as much as 15% in absolute terms as discussed in the 
text. 



Mode 


Fit Fraction (%) 


K*{892)K 


31.4 ±2.2 ±1.7 


K^(U30)K 


30.4 ±3.5 ±3.7 


K^{U30)K 


23.1 ±3.4 ±7.1 


ao{980)-K 


15.1 ±2.7 ±1.5 



of Xc -^ K*"K°TT° and that of Xc -^ K*^K^K" 
are expected to be equal. Table IIVI also shows an- 
other good agreement with the isospin expectation of 
Bixc -^ K*°K°n°)/ B{xc ^ K*°K^t:^) = 0.5 and 
B{xc -^ K*°K"tt°)/ B{xc -^ K*^TT^K^) = 0.5. 
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Figure 6: Preliminary result on 4-body decays of XcJ- 
Invariant masses of various combinations of hadrons are 
shown. 
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Table IV Branching fractions and combined error mea- 
surements for the isospin related K* Kn intermediate 
modes are listed. 



Mode 


XcO 


Xcl 


Xc2 




B.F. (%) 


B.F. (%) 


B.F. (%) 


A-*0A-±7rT 


0.56±0.15 

0.74±0.18 
0.96±0.25 


0.38±0.11 


0.59±0.14 
0.90±0.25 
0.57±0.13 
0.90±0.25 



The observation of y(4260) -^ tt^tt^ J/-0 is inconsis- 
tent with the XcJP^ molecular model 11]. Our obser- 
vation of tt^tt" J/V' rate being about half of tt+tt" J/f/; 
rate disagrees with the prediction of the baryonium 
model [12]. Evidence for the K^K~J/ip signal is not 
compatible with these two models either. Table I VII 
also shows that Y(4160) behaves very differently com- 
pared to other charmonium states above DD thresh- 
old such as ^(4040) and V'(4160) for which we set up- 
per limits in terms of cross section (cr(e+e^ —f X)) 
and branching fractions. 



3. Charmonium-like states above DD 
3.1. Y(4260) 

Y(4260) was first discovered by the BaBar Col- 
laboration via the reaction of e+e^ -^ 7^(4260) -^ 



jTr+TT-J/iP, J /ill -^ £+t 



Through the same 



production mechanism of initial state radiation, we 
also confirmed their observation based on data taken 
around the T{nS) resonances, where n is 1, 2, 3, and, 
4 Q. The invariant mass of tt+tt" J/^; based on such 
ISR production is shown in Figure [T] 
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Figure 7; Distribution of invariant mass of tt+tt^ J/tp pro- 
duced in ISR. 

We further confirmed BaBar's observation in data 
taken directly at ^ = 4260 MeV [10] at 11 cr sig- 
nificance. We also observed F(4260) -^ n'^n'^J/ip 
for the first time at significance of 5.1 a and had 
the first evidence for y(4260) -^ K+R-J/iP {3.7a). 
We measured e+e~ cross sections as a{Tr~^n~ J/ip) = 
{58t\l ± 4) pb, a{7T°Tr°J/iP) = (231^^ ± 1) pb, and 
a{K+K~J/tl)) = (9^5 ± 1) pb. Using scan data taken 
between /i of 3970 MeV and 4260 MeV, we also 
searched in 12 additional modes (transitions down to 
ip{2S), Xcj, and J/iji). No evidence of strong signals 
were seen for which we set upper limits at 90% confi- 
dence level. The results are summarized in Table IVTl 



3.2. X(3872) and Mass of neutral D 
meson 

Since X(3872) was discovered by Belle Collabora- 
tion [13| and sub seq uently confirmed by other exper- 
iments ([lJ|i[3i[l3]), many theoretical models have 
been proposed. Perhaps the most provocative theo- 
retical suggestion is that X(3872) is a loosely bound 
state of D^ and D*" mesons [17|. This idea arises 
mainly because M{D°) + M{D*°) - M{X{3872)) is 
very small. Using the average value of M(D°) of 2006 
Particle Data Group [l|, 1864.1 ±1.0 MeV, we have 
—0.9 ± 2.1 MeV for the above difference or for the 
binding energy if we assume it is a bound state of D^ 
and D*^ mesons. The large uncertainty in the differ- 
ence is partially due to the rather large uncertainty in 
mass of D^ meson. This was the motivation to mea- 
sure M{D°) more precisely using our 281 pb^^ of data 
taken at V(3770). 

We used a clean (charged particles only) mode of 
D meson decay, I?° -^ Ks(l> where Ks — * tt+tt^ 
and (j) — > K~^K~. Invariant masses of tt+tt^ and 
K^K^ are shown in Figures[8]and[9]respectively. Fig- 
ure [To] shows the invariant mass of KsK'^K" from 
which we obtained M(D°) = 1864.847 ± 0.150 ± 
0.095 MeV [H]. We then have M{D°)+M{D*°y 
M{X{3872)) = +0.6 ± 0.6 MeV. This provides a 
strong constraint for the theoretical predictions for 
the decays of X(3872) if it is a bound state of Z?° 
and D*^ mesons. The uncertainty in its binding en- 
ergy is now calling for more precise measurement on 
mass of X(3872) itself. 



4. Properties of r] 

It has been almost half a century since the rj meson 
was discovered [19]. Since then, many measurements 
have been made by many experiments. Still, almost 
all exclusive branching fractions are determined as rel- 
atives to other 77 decays. 

Based on the 27M V'(2S') sample, we measured al- 
most all the major modes (99% of generic decays of 
rj) which allowed us to determine the major branch- 
ing fractions [20|. We obtained the rj sample through 
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Table V Branching fractions (B.F.) with statistical and systematic uncertainties are shown. The symbol "x" indicates 
product of B.F.'s. The third error in each case is due to the '4){2S) — » 7Xc branching fractions. Upper limits shown 
are at 90% C.L and include all the systematic errors. The measurements of the three-hadron final states are inclusive 
branching fractions, and do not represent the amplitudes for the three-body non-resonant decays. 



Mode 


XcO 


Xcl 


Xc2 




B.F.(%) 


B.F.(%) 


B.F.(%) 




3.54±0.10±0.43±0.18 

1.48±0.13±0.18±0.08 
1.56±0.13±0.19±0.08 


1.28±0.06±0.16±0.08 

0.78±0.09±0.09±0.05 
0.78±0.09±0.09±0.05 


1.87±0.07±0.23±0.13 

1.12±0.08±0.14±0.08 
1.11±0.09±0.13±0.08 


K+K-t:"7t° 


0.59±0.05±0.08±0.03 


0.12±0.02±0.02±0.01 


0.21±0.03±0.03±0.01 


PPttV" 


0.11±0.02±0.02±0.01 


<0.05 


0.08±0.02±0.01±0.01 


K+R-'nTv'^ 


0.32±0.05±0.05±0.02 


0.12±0.03±0.02±0.01 


0.13±0.04±0.02±0.01 


K*^K^n° X K'^ -> 7r±A-o 
K'^TT^K" X if*± -. A-±7r° 


2.64±0.15±0.31±0.14 

0.37±0.09±0.04±0.02 

0.49±0.10±0.06±0.03 
0.32±0.07±0.04±0.02 
1.28±0.16±0.15±0.07 


0.92±0.09±0.11±0.06 
0.25±0.06±0.03±0.02 

0.54±0.11±0.06±0.03 


1.41±0.10±0.16±0.10 

0.39±0.07±0.05±0.03 
0.30±0.07±0.04±0.02 
0.38±0.07±0.04±0.03 
0.30±0.07±0.04±0.02 
0.42±0.11±0.05±0.03 



Table VI For each mode e^e~ -^ X, for three center-of-mass regions: the detection efficiency, e; the number of signal 
[background] events in data, A'^s [A^b]; the cross-section a{e'^e~ — > X); and the branching fraction of ?/>(4040) or i/)(4160) 
to X, B. Upper limits are at 90% CL. '-' indicates that the channel is kinematically or experimentally inaccessible. 
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38 
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58tlo±4 
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20 
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<8 


< 2 


21 
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22 
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23lg^±l 
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1 


0.07 


< 20 
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21 


3 


0.07 


9tl±l 


vJ/i' 


19 


12 


9.5 


<29 


< 7 


16 


15 


8.8 


< 34 


<8 


16 


5 


2.7 


<32 


■K^JI'4> 


23 


2 




< 10 


<2 


22 


1 




<6 


< 1 


22 


1 




< 12 


n'J/^ 








- 




11 


4 


2.5 


<23 


<5 


11 





1.5 


< 19 


7r+7r-7rOj/V> 


21 


1 




<8 


<2 


21 







<4 


< 1 


22 
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■qriJ/i) 








- 










- 
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1 




<44 


7r+7r-^(25) 








- 




12 







< 15 


<4 


19 
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■»?^(2S) 








- 
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15 







<25 


^XcO 








- 










- 




9 


11 


11.5 


<234 


7Xcl 


26 


9 


8.1 


< 50 


< 11 


26 


11 


8.7 


<45 


< 10 


26 


1 


3.3 


<30 


7Xc2 


25 


6 


8.0 


< 76 


< 17 


26 


10 


8.6 


< 79 


< 18 


27 


4 


3.3 


< 90 


TT+TT-TrOxcl 


6 







<47 


< 11 


8 







< 26 


<6 


9 







<46 


TT+TT-TT°Xc2 


4 







< 141 


<32 


8 







< 56 


< 13 


9 







<96 


7V'^7V~ (j) 


17 


26 


3.0 


< 12 


< 3 


17 


17 


6.0 


<5 


< 1 


18 


7 


5.5 


< 5 



a two-body decay oiip{2S), ^{2S) — > rjj/^ where J/^ 



+^- 



or ^^^ ). 



subsequently decays to two leptons (e+e 

That is, we have sample of about O.IM 77 decays with 

a di-lepton tag on J/4'- 

We first constrained the invariant mass of di-leptons 
to be the known mass of J/ip- We then combined 
the fitted J/ip with 77 decay products and constrained 



further to be the mass of ip{2S). In this analysis, 
the 77 decay modes we considered were 77 — > 77, Sir^, 
7r+7r^7r°, Tr^TT~^ and e+e-7. According to Ref. [l|], 



the sum of these 5 rates amounts to 99.88% of the 
total 77 decays. We then took ratios between efficiency- 
corrected yields separately for each of J/tp -^ e+e" 
and n~^iJ.~ cases in which all lepton related systematic 
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Figure 8: Invariant mass of tt+tt 
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Figure 9: Invariant mass of K'^K 



uncertainties were canceled. The resulting ratios of 77 
branching fractions are summarized in Table IVIII 

Figure [TT] shows a graphical version of comparison 
in terms of ratios of branching fractions to the single 
most precise other measurements (top of Figure [TT]) . 

By assuming that the 5 exclusive channels we con- 
sidered in this analysis cover the all r? decay modes, we 
were also able to extract absolute branching fractions 
of these 5 modes. Other possible t] decay modes are 
either now allowed and/or found to be less than 0.2% 
of generic decays 77 [J. We included 0.3% as a possi- 
ble systematic uncertainty in the absolute branching 
fraction measurements. The results are summarized 
in Table IVIIII and also shown graphically in Figure [TT] 
in terms of ratio of our branching fractions to the PDG 
2006 global fit (bottom of Figure [TTl)- Several of the 
relative and absolute branching fractions obtained in 
this analysis are either the most precise to date or first 
measurements. 

Further more, we also measured the mass of 77 me- 




M (K K+K ) (GeV) 

Figure 10; Invariant mass of KsK^ K~ . 



Table VII Ratios of ri branching fractions. For each com- 
bination, the efficiency ratio, separately for J/^j) —> e'^e~ 
and J/ip -^ /i^/i", the level of consistency between the 
J/ip — » e^e" and jJ.'^ fi~ result, expressed in units of Gaus- 



sian standard deviations, a 



! *-^^t^/ee5 



and the combined result 



for the branching ratio. The dagger symbol indicates that 
this result is most precise measurement to date. 



Channel 


eff. 
MM 


ratio 
ee 


0-,JM/ee 


branching fraction ratio 


37r'V77 


0.15 


0.15 


1.0 


0.884 ± 0.022 ± 0.019 


pi n^n /"f"f 


0.50 


0.49 


-2.2 


0.587 ± 0.011 ± 0.009t 


7r"'"7r-7/77 


0.63 


0.60 


0.2 


0.103 ± 0.004 ± 0.004+ 


e'^e~^/'y'y 


0.53 


0.52 


0.1 


0.024 ± 0.002 ± 0.001+ 


Stt^/tt+tv-t:'' 


0.30 


0.32 


2.1 


1.496 ± 0.043 ± 0.035+ 


7r+7r-7/7r+7r-7r" 


1.27 


1.24 


1.1 


0.175 ± 0.007 ± 0.006 


e e~'y/'K n^n 


1.07 


1.06 


0.5 


0.041 ± 0.003 ± 0.002+ 


e e-7/7r 7r~7 


0.84 


0.86 


0.0 


0.237 ± 0.021 ± 0.015 



son [22| . This was motivated by two recent preci- 
sion measurements that were inconsistent with each 
other. In 2002, the NA48 Collaboration reported 
M^ = 547.843±0.030±0.041 MeV [H, while in 2005, 
GEM Collaboration reported M^ = 547.311 ± 0.028 ± 
0.032 MeV 24] which was 8 standard deviations below 
NA48's resuh. 

We used the same rj sample described previously in 
this Section but used only 4 decay modes, 77 — > 77, 
37r'', TT+TT^Tr", and 7r"'"7r~7 while, again, constrain- 
ing masses of J/ip and %I}{2S). Our result, the av- 
erage of the 4 -q decay modes, is A/,, = 547.785 ± 
0.017 ± 0.057 MeV which has comparable precision 
to both NA48 and GEM results, but is consistent 
with the former and 6.5 standard deviations larger 
than the later. We note that the KLOE Collabora- 
tion also recently measured mass of the 77 meson to be 
547.873 ± 0.007 ± 0.031 MeV which was presented at 
the 2007 Lepton-Photon conference [25| . 
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Table VIII For each ri decay channel, absolute branching 
fraction measurements for J/ip —> e+e~ and J/t/j -^ M^M~ 
combined, with statistical and systematic uncertainties 
(middle column), as determined in this work. The last 
column shows the PDG fit result [3|. All but 77 are first 
measurements. 



Channel 


this work (%) 


PDG [1] (%) 


77 


38.45 ± 0.40 ± 0.36 


39.38 ± 0.26 


3^« 


34.03 ± 0.56 ± 0.49 


32.51 ± 0.28 


TT+n-n'^ 


22.60 ± 0.35 ± 0.29 


22.7 ± 0.4 


■K^-K~ 'y 


3.96 ± 0.14 ± 0.14 


4.69 ± 0.11 


e"*'e~7 


0.94 ± 0.07 ± 0.05 


0.60 ± 0.08 



5. Summary 

I have presented confirmation of BaBar's observa- 
tion of Y(4260) in di-pion transition to J/i/' along 
with a new observation through neutral di-pion tran- 
sition. Our precision measurement on M{D^) calls for 
more precise measurement on M(X(3872)). With 3M 
^(25*) sample, we have results on two-, three-, and 
four-body decays of XcJ states in which many sub- 
structures were seen in three- and four-body modes. 
Dalitz plot analyses were done for the case of 3-body 
decays. More detailed analyses can be done with the 
fuU 27M ^(25) sample. Using the 27M sample, we 
performed precision measurements on B{ri —t X) and 
M{r]). 
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Figure 11: Comparison of the results obtained in this analysis with the most precise measurements from other experi- 
ments [3, [131 (top), and the PDG 2006 global fits 1] (bottom). 



